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Abstract 

Induced  Spatial  Incoherence  (ISI)  has  been  proposed  for  KrF  laser  drivers  to  achieve  the 
high  degree  of  spatial  beam  uniformity  required  for  direct-drive  inertial  confinement  fusion. 
Although  ISI  provides  ultrasmooth  illumination  at  the  far-held  of  the  laser,  where  the  target 
is  located,  it  can  still  allow  the  beams  in  the  quasi  near-held  to  develop  time-averaged  spatial 
structure.  This  speckle,  which  arises  primarily  from  random  phase  aberration,  builds  up  as 
the  laser  beams  propagate  away  from  the  pupil  plane  located  at  the  final  amplifier  stage;  it  is 
distinct  from  any  structure  imposed  by  gain  nonuniformities  in  the  amplifiers.  Because  of  the 
spatial  incoherence,  the  speckle  is  significantly  smaller  than  that  experienced  by  coherent 
beams;  nevertheless,  it  remains  a  damage  issue,  especially  for  the  long  beam  delay  paths 
required  in  angularly-multiplexed  KrF  lasers.  This  paper  develops  a  novel  algorithm  for 
calculating  the  time-integrated  intensities,  compares  simulations  and  measurements  of  the 
near-held  speckle  in  the  Nike  KrF  laser,  and  explores  options,  such  as  aberration  reduction 
and  optical  relaying,  for  controlling  the  problem  in  future  angularly-multiplexed  KrF  drivers. 
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Introduction 


One  of  the  key  requirements  for  inertial  confinement  fusion  (ICF)  is  a  highly  uniform 
implosion  of  the  spherical  target.  In  direct-drive  ICF,  this  requires  the  spatial  profile  of 
the  focused  laser  beam  to  be  as  smooth  and  controllable  as  possible.  This  is  especially 
true  early  in  the  pulse  when  the  target  has  not  yet  developed  the  large  blowoff  plasma  that 
helps  smooth  the  spatial  structure  by  thermal  diffusion.  Irradiation  nonuniformities  at  these 
early  times  can  imprint  themselves  directly  onto  the  target,  thereby  seeding  the  Rayleigh- 
Taylor  instability  during  the  acceleration  phase.  [1] ,  [2]  A  rapid  and  effective  beam  smoothing 
technique  is  thus  required  if  the  laser  facility  is  to  be  useful  for  direct-drive  ICF. 

Induced  spatial  incoherence  (ISI)  is  a  highly  effective  technique  for  achieving  the  required 
ultrasmooth  illumination  with  a  KrF  laser  driver.  [3] -[10]  The  Nike  KrF  laser  has  produced 
smooth,  nearly  flat  ISI  focal  fluences  with  total  speckle  nonuniformities  (including  all  spa¬ 
tial  modes)  of  1%  rms  in  a  single  4  ns  beam[10]  and  as  low  as  0.15%-0.2%  when  multiple 
angularly-multiplexed  beams  are  overlapped  [8].  Using  broadband  temporally  and  spatially 
incoherent  light,  ISI  creates  the  desired  spatial  profile  at  an  object  aperture  in  the  laser  front 
end,  then  images  it  through  the  laser  system  onto  the  target  located  at  the  focal  spot.  The 
instantaneous  intensity  is  highly  nonuniform  speckle,  but  this  averages  out  to  the  desired 
profile  when  averaged  over  time  intervals  tav  much  longer  than  the  laser  coherence  time 
tc  ~  l/bandwidth.[3],[ll]  The  time-averaged  profile  remains  smooth  and  insensitive  to  laser 
phase  aberration  if  its  angular  width  is  many  times  diffraction  limit  and  large  compared  to 
the  point-spread  distribution  created  by  that  aberration. [6], [7], [8] 

Although  ISI  provides  ultrasmooth  focal  spot  illumination,  it  still  allows  significant  time- 
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averaged  speckle  structure  to  develop  in  the  quasi  near-field  of  the  laser.  This  structure, 
which  builds  up  the  beams  propagate  away  from  the  pupil-plane  images  at  the  amplifiers, 
arises  primarily  from  upstream  random  phase  aberration;  it  is  distinct  from  the  generally  long 
scalelength  variations  due  to  residual  gain  nonuniformities  across  the  amplifier  apertures.  [12] 
In  angularly-multiplexed  KrF  systems,  it  develops  primarily  in  the  long  delay  path  optics 
(e.g.,  demultiplexing  mirrors  and  focusing  lenses)  after  the  final  amplifier.  The  speckle 
nonuniformities  of  ISI  beams  are  much  smaller  than  those  created  by  coherent  light,  and  so 
far  they  have  not  caused  any  serious  damage  in  the  Nike  output  optics  because  the  average 
fluences  (typically  less  than  0.5  J/cm2)  are  well  below  damage  threshold.  Nevertheless,  they 
remain  a  potential  damage  issue  for  future  KrF  driver  designs,  where  the  pathlengths  and 
beam  sizes  would  be  similar,  but  the  average  fluences  would  have  to  be  higher.  [13]  A  fast 
and  reliable  algorithm  is  therefore  needed  to  simulate  and  minimize  the  near-field  speckle 
problem  in  these  future  designs.  In  this  paper,  we  derive  such  an  algorithm,  compare  its 
simulations  with  measurements  on  the  Nike  laser,  and  use  it  to  explore  options  for  controlling 
the  problem  in  future  angularly-multiplexed  systems. 

We  examine  two  techniques  for  calculating  the  average  intensities.  The  first  models  the 
instantaneous  stochastic  optical  amplitude  at  the  object-plane  by  a  2D  Gaussian-distributed 
random  complex  array,  then  propagates  it  through  the  laser  system  (whose  complex  final 
aperture  transmission  includes  the  phase  aberration)  to  evaluate  the  instantaneous  intensity. 
One  then  repeats  this  procedure  with  multiple  new  statistical  realizations  to  accumulate 
the  average  intensity.  The  second  technique,  which  is  much  faster,  directly  evaluates  the 
average  intensity  envelope  at  both  the  near-field  or  far-field  by  using  a  generalized  version  of 
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partially-coherent  imaging  theory. [11], [14], [15]  Starting  with  the  incoherent  object  envelope, 
it  calculates  the  optical  autocorrelation  function  before  the  final  amplifier  aperture  from 
the  Van  Cittert-Zernike  theorem,  then  combines  it  with  the  phase-aberrated  transmission 
to  derive  convolution  integrals  that  relate  the  downstream  intensity  envelopes  to  the  object. 
(For  the  near-held,  the  point-spread  distribution  is  an  aberrated  Fresnel  diffraction  formula.) 
The  convolutions  are  then  evaluated  by  Fast  Fourier  Transform  (FFT)  methods.  Benchmark 
tests  show  excellent  agreement  (within  <  1%  rms)  between  the  two  techniques. 

We  then  compare  the  theory  with  calibrated  photochromic  film  images  of  two  Nike  laser 
beams.  The  simulations  assume  uniform  illumination  at  the  60  cm  final  amplifier  aperture, 
a  uniformly  illuminated  circular  object  aperture  of  width  75  x  Diffraction  Limit  (75  XDL), 
and  random  phase  aberration  of  ~  A/4  rms  amplitude  and  ~  15  XDL  far-field  point-spread 
width.  In  good  qualitative  agreement  with  the  measurements,  the  simulations  show  that 
as  the  distance  from  the  amplifier  increases,  (a)  the  structure  becomes  coarser  and  (b)  the 
hot-spot  amplitudes  increase  at  first,  then  appear  to  saturate  and  slowly  decrease.  In  all 
cases,  the  measured  and  simulated  far-field  profiles  remained  smooth.  [10] 

Finally,  we  perform  additional  simulations  to  explore  options  for  reducing  the  near-field 
speckle.  One  possibility  is  to  reduce  the  phase  aberration;  e.g.,  a  two- fold  reduction  in 
point-spread  width  also  reduces  the  intensity  variations  by  ~  two-fold.  The  second  option 
would  reduce  the  hot-spots  by  using  optical  relaying,  but  that  would  probably  be  expensive 
to  implement  in  some  of  the  longer  beam  paths  required  for  demultiplexing. 

Although  this  paper  addresses  only  the  case  of  ISI  in  angularly-multiplexed  KrF  lasers, 
the  problem  of  quasi  near-field  structure  is  not  limited  to  such  systems.  Other  beam- 
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smoothing  techniques,  such  as  smoothing  by  spectral  dispersion  (SSD) ,  [16] ,  [17]  optical  fiber 
smoothing  (OFS),  [18]  and  partially  coherent  light  (PCL),[19]  are  also  subject  to  time- 
integrated  near- held  structure  due  to  laser  phase  aberration.  These  techniques  are  normally 
applied  to  glass  lasers,  where  optical  relaying  is  easier  to  implement  because  there  are  fewer 
beams;  however  glass  lasers  are  also  more  susceptible  to  small-scale  self-focusing,  which  tends 
to  amplify  the  near-held  nonuniformities. 

Theory 

Figure  1  illustrates  the  implementation  of  ISI  in  the  Nike  laser.  A  spatially  and  tem¬ 
porally  incoherent  oscillator  uniformly  illuminates  the  object  aperture,  whose  transmission 
produces  the  desired  time-averaged  transverse  spatial  profile  ( Iq  (xo))  =  (| Eq  (xo,£)|  )  from 
the  rapidly- varying  complex  speckle  amplitudes  Eo  (xo,  t ).  This  light  is  relayed  through  the 
laser  system,  then  focused  to  produce  the  image  (/p(x^))  =  (\EF  (x^,  t)|2)  at  the  target. 
Each  amplifier  is  placed  at  or  near  an  image  of  the  object’s  Fourier-transform  (FT)  plane, 
where  the  light  from  each  point  of  (Jo  (xo))  encodes  itself  into  a  bundle  of  rays  that  traverse 
the  amplifier  at  the  same  unique  held  angle.  Because  these  angles  are  small,  each  bundle  is 
amplified  the  same  and  (IF  (x^))  remains  insensitive  to  laser  phase  and  gain  nonuniformities. 
This  insensitivity  can  also  be  explained  by  an  alternative  picture  in  which  the  object  prohle 
encodes  itself  into  Np  iT  1  small  coherence  zones  in  the  FT  plane;  [4]  as  long  as  these  zones 
remain  small  compared  to  the  nonuniformity  scalelengths,  they  are  not  strongly  affected. 

For  KrF  bandwidths  Azz  ~  1  THz  and  a  carrier  frequency  u0  =  c/A  ~  1200  THz,  we 
can  usually  treat  the  light  propagation  through  the  laser  system  as  quasi-mono  chromatic. 
This  approximation  is  justihed  as  long  as  the  time-averaged  near-held  structure  arises  from 
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interfering  rays  whose  optical  path  differences  are  less  than  the  optical  coherence  length 
ctc  ~  c/Au  ~  300  pm.  In  the  Appendix,  we  show  that  this  condition  is  well  satisfied  in 
nearly  all  the  simulations  presented  here  because  the  small  transverse  width  of  the  spatial 
coherence  zones  effectively  limits  the  path  differences;  the  only  marginal  case  is  that  of 
spatially  coherent  (but  still  finite  Av)  light.  We  will  thus  simplify  the  notation  by  dropping 
explicit  references  to  the  time  variation  of  Eo  (xo,t)  and  all  other  complex  amplitudes. 

In  the  Nike  laser,  the  object’s  FT  plane  is  located  at  a  rectangular  entrance  pupil  aperture 
of  width  DE  located  a  focal  distance  /#  beyond  the  input  lens.  The  lens  collimates  some  of 
the  light  from  the  incoherent  object  to  overfill  this  aperture,  which  is  then  imaged  into  the 
amplifiers.  The  instantaneous  speckle  amplitude  is  thus [20] 

Ee  (xb)  =  (XfE)-1  d  [Eo,  -xE/  (A fE)\  ,  (1) 

where  A  is  the  optical  wavelength  and 

3  (E,  ±x)  =  J  E  (£)  exp  (±i27r£  •  x)  d2£  (2) 

defines  the  2D  spatial  FT  of  E  (£).  Because  of  the  overfill,  we  can  assume  approximately 
uniform  average  intensity  (\Ee  (x^)]  )  within  the  aperture. 

The  laser  system  splits  Ee  (x®)  into  numerous  angularly- multiplexed  beams,  then  magni¬ 
fies  and  image-relays  each  one  through  a  chain  of  amplifiers  to  the  final  stage,  whose  aperture 
Da  acts  as  the  optical  stop  or  exit  pupil.  In  the  absence  of  laser  gain  and  phase  aberration, 
the  amplitude  just  before  this  aperture  would  then  be  Ea  (x^)  =  Ee  (x^/Mg^)  / MEa ,  where 
ME a  is  the  geometrical  optics  magnification  MEa  (A  Da/De  >  1,  and  the  average  intensity 
(\Ea  (Xj4)|2)  is  again  uniform.  Because  of  the  amplifier  imaging,  one  can  approximately  lump 
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all  the  laser  gain  and  phase  nonuniformities  into  a  single  complex  aperture  transmission 


Ta  (xa)  =  A  (xA)  y/ G  M  exp  (xA)]  ,  (3) 

where  G*  (xA)  is  the  total  power  gain,  (j)  (xA)  is  the  total  phase  aberration,  and  A  (xA)  is  the 
hard  aperture  transmission.  For  a  square  aperture,  A  (xA)  =  1  when  \xA\  ,  \yA\  <  DA/2  and 
0  otherwise.  The  amplitude  just  beyond  the  aperture  is  thus 

E'a  (xa)  =  Ta  (xa)  Ea  (xa)  =  Ta  (xa)  Ee  (xa/Mea)  /Mea  .  (4) 

The  spatial  variation  of  G  (xA)  includes  the  effect  of  gain  saturation  as  well  as  laser  excitation 
nonuniformity,  and  is  assumed  to  be  only  long  scalelength.  In  principle,  rapidly-varying 
speckle  could  create  short  scalelength  gain  nonuniformities  due  to  saturation,  but  the  changes 
occur  so  rapidly  that  the  KrF*  population  is  unable  to  respond.  [4]  Saturation  ensures  that 
the  average  intensity  (|i?A(xA)|2)  is  similar  to  the  final  amplifier  excitation  profile,  which 
would  be  approximately  uniform  in  a  practical  system.  The  hard  aperture  diffraction  created 
by  A(xa)  will  cause  a  slight  blurring  of  the  target  beam  profile  (Jp  (xf)),  but  the  effect  is 
negligible  as  long  as  the  width  of  ( IE  (xp))  remains  many  XDL.  In  the  Nike  laser,  this  profile 
is  typically  ~  75  XDL;  in  a  fusion  driver  it  is  likely  to  be  even  larger. 

One  can  write  the  complex  amplitudes  in  a  compact  form  suitable  for  numerical  analysis 
by  rescaling  transverse  near-held  coordinates  to  the  geometrical  optics  beamwidths,  the 
far-held  (object  and  image)  coordinates  to  the  diffraction-limited  widths,  and  propagation 
distances  to  Fresnel  lengths;  e.g.,  xA  scales  to  DA ,  while  xq  scales  to  do  =  MEAfE\/DA : 

XA  =  xA/ Da  ,  X0  =  xq/ do  . 
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Expressions  (l)-(4)  then  give  the  amplitude  E'A  (X^)  =  Ta  (X^)  Ea  (X^) ,  where  the  complex 
aperture  transmission  now  satisfies  Ta  (X^)  =  0  for  \Xa  \  ,  |Ya|  > 

Ea  (Xa)  =  Fa  J  Eo  (X0)  exp  (— ®2tj tXa  •  X0)  d2X0  =  FA  $  (E0,  -XA)  ,  (5) 

and  Fa  =  MeaJe^/ Da  <C  1  is  an  inverse  Fresnel  number.  With  these  scalings,  all  near-held 
beamwidths  are  limited  to  ~  1,  but  the  far-held  widths  can  range  over  many  XDL. 

Figure  2  shows  one  of  the  angularly-multiplexed  beams  at  the  output  of  the  optical 
system.  Just  after  the  common  aperture  stop  A,  a  weak  focusing  optic  converges  each 
beam  to  its  collimator  mirror  C  at  distance  zAc-  The  collimated  beam  of  geometrical  optics 
width  Dq  <  Da  then  propagates  over  a  delay  distance  zcl  to  the  final  focusing  lens  L, 
which  is  still  in  the  quasi  near-held  of  A.  The  length  of  the  delay  path  is  controlled  by 
the  location  of  demultiplexing  mirror  M,  and  is  different  for  each  beam.  In  the  case  of 
interest,  where  the  beam  underfills  all  the  optical  elements  located  downstream  from  A, 
we  can  calculate  its  propagation  directly  from  E'A  (XA,t).  We  have  assumed  that  the  phase 
aberration  (j)  {XA)  in  Ta  (X^)  adequately  represents  all  of  the  contributions  introduced  by  the 
entire  laser  optical  system,  including  those  from  the  collimation  and  time  delay  optics.  This 
approximation  is  reasonable  because  most  of  the  aberration  originates  in  or  around  the  large 
aperture  amplifiers,  where  the  small  held  angles  associated  with  the  large  beamwidth  are 
more  susceptible  to  perturbations  by  the  optics  and  air  paths.  Phase  perturbations  beyond 
the  collimators  will  be  less  effective  in  producing  speckle  at  the  focusing  lenses  because 
the  held  angles  are  larger  and  the  remaining  propagation  distances  are  smaller;  thus  the 
simulations  presented  here  will  slightly  overestimate  the  structure. 

We  calculate  the  beam  propagation  using  the  Huygens-Fresnel  integral.  [20]  In  the  scaled 
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coordinates  Xa  =  and  Xc  =  x.c/Dc ,  the  converging  phase  curvature  imposed  by 

the  weak  focusing  optic  at  A  is  —  n  |.DaXa|2  /  (A/a)  ,  while  the  diverging  curvature  imposed 
by  the  collimator  C  is  +7T  l-D^Xd2  /  (A  \fc\)  •  The  focal  lengths  /a  and  fc  are  related  to  the 
propagation  distance  zAc  and  geometrical  optics  beamwidths  by 

\fc\/fA  =  Dc/DA,  \fc\  +  zAc  =  /a  • 

In  the  paraxial  approximation,  the  collimated  beam  amplitude  just  beyond  C  is  then 

„  ,  D\  j .  |DcXc|2j  /•„  ,  (  .  |^XA|2j 

Bc(Xc)  =  {-'*—) 
x  exp  f 1-DcXc  -  .DaXaI2)  d2XA 

V  A^AC  / 

=  M  V  [  E'a  (XA)exp  (W-  |XC  -  X^l2)  <i2X4,  (6) 

Mac ^ ac  J  \^ac  J 

where  Mac  =  Dc/Da  =  |/c|//a  is  the  spatial  magnification  (now  <  1)  and  ZAc  = 
A zac/  ( DADc )  1  is  the  propagation  distance  zac  scaled  to  the  Fresnel  length  DADc/ A. 
Similarly,  the  amplitude  just  before  the  final  focusing  lens  L  at  distance  zql  is  given  by 

El  (Xi)  =  J^lJEc  (Xc)  exp  (£ 7 |Xi "  Xc|2)  d2Xc  ’ 

where  XL  =  xL/F>L  =  xL/Dc  and  ZCL  =  AzCl/  (DcDl)  =  \zCl/Dc  (Mcl  -»•  1)  •  We  can 
also  relate  El  directly  to  E'A  by  combining  these  expressions  and  integrating  over  X,^ : 

El  (Xt)  =  1  f  E'a  (Xa)  exp  (/^  |XL  -  X4|2  j  d2X4 ,  (7) 

Mac^al  J  \^al  J 

where  Zal  =  ZAc  +  -Zol  .  Because  these  expressions  have  the  same  form,  one  can  describe 
all  the  quasi  near-held  propagation  simply  by  summing  the  appropriate  dimensionless  path- 
lengths  Zac  ,  Zqli  etc,  even  if  the  beam  expands  or  converges  along  the  path. 
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The  final  lens  provides  the  phase  curvature  — 7T  \  DlX.l\  /  (A f£)  required  to  focus  El  to 
its  far-field  amplitude  EF.  Substituting  E'L  (XL)  =  EL  (XL)  exp  [—in  | id>^X ^  | 2  /  (A Jl)\  into 
the  Huygens-Fresnel  integral  for  a  distance  equal  to  the  focal  length  ,  we  obtain 

EF  (Xf)  =  (1  /Fl)  exp  (in tFl  |Xf|2)  £  [EL,  -XF]  ,  (8) 

where  X^  =  x.F/dF  is  scaled  to  the  diffraction-limited  width  dF  =  JlX/Dl  and  Fl  = 
X/l/D^.  As  in  the  near-field  case,  we  can  write  EF  directly  in  terms  of  E'A  by  substituting 
(7)  and  integrating  over  XL: 

Ef  (Xf)  =  (1  /MacFl)  exp  [nr  (FL  -  ZAL)  |XF|2]  3  (E'A,  -XF)  .  (9) 

This  recovers  the  well-known  result  that  the  far-field  amplitude  differs  from  the  exact  Fourier 
transform  of  the  aperture  amplitude  only  by  a  phase  curvature  term,  which  does  not  affect 
the  focal  intensity  distribution  IF  (XF,t)  =  \EF  (X^)|2  . 

The  above  expressions  can  be  evaluated  numerically  by  FFT  methods.  For  example, 
FFT  calculates  the  far-field  EF  (XF)  directly  from  Eqs.  (8)  or  (9).  In  principle  one  can  also 
use  FFT  to  evaluate  El  (Xl)  directly  by  expressing  the  integral  in  (7)  as  $  (EA,  —Xl/Zal)  , 
where  EA  (X^)  =  E'A  (XA)  exP  [}^  ( | X^4 1 2  +  |XL|  )  /ZAl\  ;  however  that  approach  is  not  well 
suited  to  the  near-field  regime  because  the  small  values  of  ZAl  introduce  rapid  phase  os¬ 
cillations  in  EA(XA).  Instead,  we  use  (7)  to  relate  the  Fourier-transformed  amplitudes 
El  (f3)  =  $  (El,  ~P)  and  E'A  (/3)  =  $  (EA,  — /3)  by  a  unitary  propagator;  i.e., 

El  (P)  =  (Mac)~1  E'a  (p)  exp  (  inZAL  \P\2)  (10) 

The  inverse  transform  El(Xl)  =  $  ^El,+X.lJ  then  yields  the  desired  amplitude,  which 
gives  the  speckle  intensity  IL  (XL,t)  =  \EL  (Xi)|2. 
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The  time-averaged  ISI  propagation  can  always  be  simulated  by  superimposing  the  random 
speckle  patterns  calculated  from  a  large  ensemble  of  statistically-independent  realizations  of 
the  complex  object  amplitude  Eq  (x).  In  each  realization,  Eq  (x)  is  modelled  by  a  2D  circular 
array  of  random  complex  numbers  satisfying  Gaussian  statistics  with  equal  RMS  values.  This 
produces  a  ”  top-hat”  average  intensity  distribution  imbedded  in  an  otherwise  empty  square 
window  to  provide  guard  bands.  For  each  realization,  one  would  use  the  FFT  to  calculate 
the  exit  pupil  amplitude  EA(XA)  from  Eqs.  (4)  and  (5),  assuming  a  uniform  intensity 
transmission  T  (X^)|  =  1  in  a  square  aperture  and  a  random  phase  aberration  (j)  (X^)  of 
specified  RMS  amplitude  and  gradient.  (This  aperture  is  also  imbedded  in  a  larger  window  to 
provide  guard  bands.)  To  allow  for  overfilling  at  the  aperture,  one  must  choose  the  number  of 
points  across  the  object  to  be  somewhat  larger  than  the  width  (in  XDL  units)  of  the  desired 
beam  profile  at  the  target.  The  far-held  and  quasi  near-held  intensity  prohles  are  then 
evaluated  using  the  formulas  derived  above.  Although  this  approach  is  straightforward,  the 
large  number  of  realizations  required  to  adequately  smooth  out  the  statistical  nonuniformities 
leads  to  excessive  calculation  times,  even  for  moderately  structured  beams. 

In  the  simulations  presented  here,  we  use  a  faster  alternative  approach  that  relates  the 
average  intensity  prohle  (i.e.,  the  intensity  envelope)  directly  to  the  spatial  autocorrelation 
function  just  after  the  exit  pupil  at  A: 

T\  (Xa,  X’a)  =  (E’A  (Xa)  E’X  (XI a))  =  Ta  (Xa)  TJ  (X'a)  Ta  (xa,  x'a)  .  (11) 

where  ()  denotes  an  ensemble-average.  The  incident  correlation  function  T^X^X^)  = 
(EA  (X^)  E*a  (X!a))  is  related  to  the  correlation  function  To  (Xo,X'0)  =  (Eq  (Xq)  Eq  (X'0)) 
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at  the  object  plane  by  expression  (5): 


rA(xA,x'A)  =  F* 


r0  (Xo,  x'0)  exp  [-fix  (X,  •  Xo  -  x'A  ■  x;,)]  <f2x0<;Jx' 


O  ■ 


Because  average  intensity  nonuniformities  are  lumped  into  \TA  (X^)|2  instead  of  (X^,  X^) , 
we  can  write  (X^,  X'A)  ~  (X^  —  X'A).  One  can  model  this  condition  and  still  allow  a 

nonuniform  object  intensity  envelope  (Jo  (Xo))  by  adopting  the  incoherent  approximation 
To  (Xq,  X'o)  ~  (Jo  (Xo))  <5  (Xq  —  X'G)  ,  thereby  giving  the  Van  Cittert-Zernike  result  [11] 


VA  (XA  -  X'A)  =  Fj  5  [(/o>  ,-(XA-  X(4)J  . 


(12) 


We  now  construct  a  numerical  algorithm  based  on  FFT  techniques  that  directly  evaluates 
the  envelope  profile  (I  (X))  =  <(|£'(X)|2)  at  either  the  focal  spot  [(/  (X))  — >  (IF  (Xi?))]  or 
the  quasi  near- held  [e.g.,  (I  (X))  — >  (IF  (X^))]  by  modifying  a  well-known  formalism  used 
to  describe  partially-coherent  images. [11], [14]  Combining  expression  (11)  with  either  (9)  or 
(7),  then  substituting  (12)  we  can  write  the  general  result 

(I  (X))  =  f J  rA  (XA,  X'A)  Ta  (Xa)  K  (X,  XA)  T\  (X'J  K*  (X,  X'a)  d2XAd2X'A 


J  (Jo(Xo))  Fa  J  TA(XA)K(X,XA)exp(-i2nXA-X0)d2X/ 


d2X0,  (13) 


where  the  respective  far- held  and  near- held  kernels  K  — >  Kf,l  are: 


KF(XF,XA)  = 

( MacFl )  1  exp  (-i2nXF  ■  XA) 

(14) 

Kl  (Xl,  Xa)  =  (MacZal)-1  exp  (mZA\  |XL  -  XA|2)  . 

With  these  kernels,  expression  (13)  reduces  to 

(15) 

(If(Xf))  =  j 

f  (Io  (X0))  SF  (XF  +  X0)  d2X0 

(16) 

(II  (Xl))  =  j 

f  (Io  (Xo))  SL  (XL  +  ZALXD )  d2X0  , 

(17) 
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where  the  point-spread  functions  are 

Sf (X)  =  \(C/Fl)  5(T4,-X)|2  (18) 

Si(X)  =  f  TA(XA)eKP(^-\X-XA\2)d2XA  (19) 

^AL  J  \^AL  J 

and  C  =  Fa/ Mac  ■  Equation  (16)  is  the  well-known  convolution  that  relates  the  focal  image 
profile  to  the  incoherent  object  and  a  point-spread  function  Sp  (X)  dominated  by  the  phase 
aberration  in  Ta  (Xyi).  As  long  as  Sf  (X)  remains  narrow  compared  to  (Jo  (Xo)),  its  detailed 
shape  and  structure  will  have  little  effect  on  (Ip  (X.p)),  except  for  some  rounding  near  the 
edges.  This  is  illustrated  in  Fig.  3,  which  simulates  ISI  on  the  Nike  laser  under  typical 
operating  conditions.  The  near-held  result  (17)  is  equivalent  to  the  convolution  used  to 
describe  Fresnel  diffractive  effects  in  radiography,  [15]  but  its  point-spread  function  (19)  has 
now  been  generalized  to  include  phase  aberration.  Expression  (17)  is  dominated  by  the 
structure  in  Sl  (X)  when  Z4Lmax(Xo)  <C  max(Xx);  it  begins  to  smooth  out  only  when 
Zap  max  (Xq)  becomes  comparable  to  max  (XL) .  (As  expected,  (17)  approaches  the  far-held 
shape  (Ip)  in  the  Fraunhofer  limit  Zal  —■ "  oo.) 

The  integrals  in  (18)  and  19  are  identical  to  those  of  the  respective  far-held  and  near-held 

amplitudes  (9)  and  (7),  except  for  replacing  E'A  (X^)  by  Ta  (X^) .  Their  evaluation  therefore 

follows  the  same  steps  as  those  for  the  amplitudes;  i.e.,  Sp  (X)  is  found  by  applying  the  FFT 

directly  to  expression  (18),  but  Sp  (X)  requires  that  we  express  Ta(X-a )  hr  terms  of  its 

Fourier  transform  Ta  (/3)  =  $  {Ta,  —(3)  before  integrating  (19)  to  obtain 

/2  2 
Tl  {(3)  exp  (*27tX  ■  (3)  d2[3  =|c7ff(fL,x)|  ,  (20) 

where  TL{(3 )  =  Ta  {(3)  exp  (—ittZal  \(3\ 2)  is  calculated  from  the  same  unitary  propagator 
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used  for  the  complex  amplitudes  in  Eq.  (10). 

Expressions  (16)  and  17  can  be  efficiently  solved  numerically  by  taking  the  Fourier  trans¬ 
forms  (^1  (/3)\  =  3  ((/) ,  —(3),  applying  the  convolution  theorem[ll],  and  recalling  (12): 

=  3((Io),+(3)3(Sf,-(3)  =  FX2Ta(-[3)  3(Sf,-P)  (21) 

(h(P))  =  3({Io),+Zal(3)3(Sl,-(3)  =  FX2Ta(-Zal(3)  $(Sl,-0)  (22) 

The  envelopes  then  follow  from  the  inverse  transforms  (If,l  (X-f,l))  =  3  (Jf,l  ,  +Xpj^ . 

Figure  4  compares  the  envelope  prohles  with  the  intensities  averaged  over  8000  statisti¬ 
cal  realizations  at  both  the  near-  and  far-held  planes.  Even  with  the  reduced  widths  and 
structure  in  these  simulations  (e.g.,  32  XDL  instead  of  75  XDL),  the  averaging  approach 
required  several  hours,  whereas  the  direct  envelope  calculations  took  less  than  30  seconds. 
This  large  difference  in  running  times  and  the  excellent  agreement  at  all  locations  confirm  the 
value  of  the  new  algorithm.  Except  for  the  residual  Fresnel  diffraction  spikes  at  the  edges, 
the  envelope  structure  at  the  focusing  lens  arises  entirely  from  phase  aberration;  when  this 
aberration  is  switched  off,  only  the  Fresnel  spikes  remain. 

The  above  analysis  was  carried  out  for  the  Nike  baseline  case  shown  in  Fig.  2,  where  a 
diverging  mirror  collimates  the  beam  before  it  focuses,  but  it  can  apply  equally  well  to  the 
alternative  case  where  a  converging  mirror  collimates  the  beam  after  it  traverses  a  focus. 
(See  Fig.  10)  This  option  would  require  the  output  beams  to  be  enclosed  in  an  evacuated 
tunnel,  which  would  probably  be  desirable  in  any  case  to  avoid  excessive  phase  distortion 
from  air  turbulence  or  nonlinear  optical  effects  in  the  long  demultiplexing  paths.  To  treat 
this  configuration,  the  above  formalism  requires  only  two  modifications:  (a)  change  the 
sign  of  the  scaled  propagation  distance  ZAc ,  and  (b)  invert  all  prohles  at  and  beyond  the 


14 


collimators  [e.g.,  (Ic  (XL))  — >•  (Ic  (-XL))  . 


Results  and  Discussion 

In  this  section,  we  compare  ISI  average  intensity  simulations  with  measured  quasi  near¬ 
field  speckle  on  the  Nike  laser,  compare  the  ISI  speckle  with  that  of  spatially-coherent  light, 
and  explore  options  for  reducing  the  structure.  (Detailed  comparisons  of  Nike  far-field  profiles 
are  described  in  refs.  [7]  and  [10].)  The  simulations  use  the  same  parameters  as  in  Fig.  3; 
i.e.,  a  ”  top-hat”  object  profile  of  diameter  NF  =  75  XDL,  a  uniformly  illuminated  square 
amplifier  aperture  of  width  Da  =  60  cm,  and  random-phase  aberration  of  A/4  rms  amplitude 
characterized  by  a  point-spread  function  of  ~  15  XDL  half-energy  width.  In  all  cases,  the 
collimated  beamwidths  are  D(:  =  15  cm  and  the  aperture-collimator  distances  are  Zac  =  53 
m,  thus  giving  the  scaled  A-C  distance  Zac  =  A zac /  ( DaDc )  =  1.46  x  10~4. 

Using  a  setup  similar  to  that  described  in  ref.  [12],  we  recorded  time-integrated  near-field 
images  of  several  Nike  beams  on  photochromic  UV-sensitive  PERM  film  (manufactured  by 
GAF  Chemicals). [21]  For  each  beam,  we  measured  the  image  at  both  the  collimator  and 
the  target  chamber  window  during  the  same  laser  shot  by  placing  film  over  one  half  of  the 
recollimator  and  allowing  the  unblocked  half  of  the  beam  to  propagate  on  to  the  window. 
(See  Fig.  2.)  The  windows  are  located  one  meter  beyond  the  6  meter  focusing  lenses,  so  the 
geometrical  beamwidths  at  the  windows  are  Dw  =  5/6 Dc  =  12.5  cm.  To  avoid  saturating 
the  film,  we  fired  the  system  only  up  to  the  20  cm  driver  amplifier,  but  not  the  60  cm 
final  stage.  This  does  not  affect  the  measurements  because  the  quasi  near-field  structure 
is  dominated  by  imperfections  in  the  final  amplifier  optics,  and  does  not  depend  on  laser 
intensity.  The  film  was  digitized  using  a  commercial  8-bit  scanner  (HP  ScanJet  lie).  We 
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calibrated  the  film  response  at  248  nm  by  using  repeated  shots  from  a  small  KrF  laser  to 
vary  the  total  film  exposure  from  1  to  50  mJ/cm2,  thus  creating  a  custom  greyscale  target. 
A  scan  of  this  target  revealed  that  its  weak  saturation  behavior  was  well  approximated  by 
an  exponential  function,  which  was  then  used  to  linearize  the  raw  image  data. 

Figure  5  shows  the  recorded  ISI  images  and  lineouts  at  the  recollimator  mirror  array 
and  target  chamber  windows,  while  Fig.  6  shows  the  corresponding  simulations.  (In  these 
and  in  all  other  near-field  lineouts,  the  fluences  are  normalized  to  the  average  value  across 
the  geometrical  optics  beamwidth  between  —0.5  and  +0.5.)  The  measured  fluences  in  Fig. 
5  lie  well  within  the  linear  portion  of  the  PERM  film  response,  with  average  intensities 
~  5  mJ/cm2.  Beam  #38  (Figs.  5b  and  6b)  has  one  of  the  shorter  collimated  propagation 
paths,  with  a  total  collimator- focusing  lens  distance  zcl  =  41  meters  and  a  scaled  distance 
ZCl  =  A zcl/ D2C  =  4.52  x  10-4;  beam  #10  (Figs.  5c  and  6c)  has  one  of  the  longer  paths, 
with  a  total  distance  zcl  =  76  meters  and  Zcl  =  ^zcl/D2c  =  8.38  x  10  4.  Because  the 
aperture-recollimator  distance  is  nearly  the  same  for  all  beams,  we  show  the  recollimator 
results  only  for  beam  #10.  Figures  5  and  6  clearly  show  that  the  structure  becomes  coarser 
as  the  light  propagates  farther  from  the  amplifier  aperture,  as  expected  from  Eq.  (17). 
Figure  6  also  indicates  that  the  nonuniformity  amplitudes  increase  at  first,  then  saturate 
and  slowly  decrease  as  the  pathlength  increases. 

Figure  7  compares  the  RMS  transverse  mode  amplitudes  corresponding  to  the  structure 
shown  in  Figs.  5  and  6.  To  obtain  these  spectra,  we  used  a  procedure  similar  to  that  of  ref. 
[10];  i.e.,  we  averaged  the  ID  FFT  power  spectra  over  a  large  number  of  horizontal  lineouts, 
then  took  the  square  root.  For  the  simulated  data,  the  lineouts  included  the  entire  beams; 
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for  the  measured  data,  we  averaged  over  150  lines,  excluding  those  in  the  weaker  portions 
of  the  beams  and  those  containing  the  anomalously  large  hot  spot  at  the  top  center  of  the 
image  in  Fig.  5b.  To  avoid  spurious  contributions  due  to  rapid  rolloff  at  the  edges  of  the 
beams,  we  applied  a  ID  Hanning  window [22]  to  the  data.  (This  affects  only  the  first  nonzero 
mode.)  The  RMS  amplitudes  are  scaled  to  the  zeroth  component,  and  the  mode  numbers 
correspond  to  the  number  of  cycles  across  the  geometrical  optics  beamwidth. 

Figures  7b, c  show  good  agreement  between  the  measured  and  simulated  data  except 
for  the  broadband  shelf  above  20  cycles,  which  appears  to  be  him  noise.  A  comparison 
between  Figs.  7b  and  7c  also  shows  the  spectral  narrowing  due  to  the  suppression  of  higher 
spatial  frequencies  in  the  longer  propagation  path  of  beam  #10,  as  noted  above.  Figure 
7a  shows  generally  lower  amplitudes  at  the  recollimators,  as  clearly  seen  in  Fig.  6,  but  the 
measured  beam  profile  contains  additional  higher  spatial  frequency  structure  not  present  in 
the  simulated  profile.  This  structure,  which  also  shows  up  as  the  weak  semicircular  rings 
and  radial  striations  in  Fig.  5a,  appears  to  be  caused  by  HR  optical  coating  imperfections 
in  a  large  circular  mirror  located  between  the  20cm  and  60cm  amplifiers.  Because  it  occurs 
predominantly  at  spatial  frequencies  above  those  in  the  simulated  image,  it  tends  to  diffract 
away  more  effectively  with  increasing  propagation  distance,  and  thus  has  little  impact  on 
the  results  shown  in  Figs.  7b, c. 

Figure  8  shows  the  simulated  structure  in  beam  #10  for  the  case  of  spatially-coherent 
light.  This  corresponds  to  a  point  source  object  ( Io  (Xo))  — >  6  (Xo),  so  (#  (X#)  in  Eq. 
(17)  reduces  to  the  usual  Fresnel  diffraction  result  (19).  Comparing  Fig.  8  with  the  ISI  case 
in  Figs.  6a  and  6c,  we  see  that  the  structure  is  nearly  the  same  at  the  recollimator,  with  only 
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a  small  enhancement  of  the  Fresnel  spikes  at  the  edges.  As  the  propagation  path  increases, 
however,  the  structure  becomes  far  more  severe  than  that  of  the  ISI  case,  with  hot-spots  as 
large  as  10  times  the  average  fluence.  The  difference  occurs  because  ISI  filtered  out  more 
of  the  high  spatial  frequency  components  with  increasing  propagation  path.  This  can  be 
easily  explained  by  using  either  Eq.  (17)  or  the  alternative  ISI  model  mentioned  above,  in 
which  the  object  encodes  itself  into  small  independent  coherence  zones  of  width  cIa  =  Da/Nf 
at  the  amplifier  aperture.  (NF  =  75  XDL  in  these  simulations)  Because  ISI  averages  out 
interference  from  aperture  points  separated  by  more  than  (1a  ,  it  effectively  cuts  off  all  of  the 
structure  with  transverse  wavelengths  shorter  than  some  minimum  value  A.  To  estimate  A 
at  the  collimator,  we  note  that  according  to  Eq.  (6),  any  two  aperture  points  and  X'A 
separated  by  scaled  distance  (1a/ Da  will  produce  an  interference  pattern  of  scaled  minimum 
wavelength  A c/Dc  =  A TfZac]  a  similar  argument  applied  to  Eq.  (7)  gives  the  analogous 
result  Af/Df  =  NfZal  =  AfF  ( Zac  +  Zcl)  at  the  lens.  Recalling  the  definitions  of  Zac  and 
Zee-,  we  then  obtain  Ac/ Dc  —  1/91.3  and  A F/DF  ~  1/13.6,  which  is  consistent  with  Fig. 
8.  At  the  collimator,  ISI  allows  components  with  up  to  91  cycles  across  the  beamwidth,  so 
it  reduces  only  the  sharp  Fresnel  spikes  at  the  edges;  at  the  lens,  however,  ISI  allows  fewer 
than  14  cycles,  so  it  suppresses  almost  all  of  the  sharp  high  spatial  frequency  structure.  The 
ISI  filtering  process  is  less  effective  in  beam  A 38.  allowing  up  to  22  cycles  across  the  beam 
at  the  lens  plane  (A F/DF  ~  1/22.3).  This  contrast  is  also  evident  in  Fig.  7. 

One  obvious  way  of  reducing  the  quasi  near-held  structure  is  to  reduce  the  phase  aberra¬ 
tion,  as  illustrated  in  Fig.  9.  Comparing  Figs.  9  and  6a, b,  we  see  that  a  ~  2 -fold  reduction 
in  point-spread  width  would  also  reduce  the  nonuniformities  by  ~  2-fold.  A  reduction  of 
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this  magnitude  should  be  possible  in  future  KrF  drivers,  which  would  use  evacuated  beam 
paths  and  take  advantage  of  ion  milling  techniques  to  improve  optical  surface  quality.  It  is 
also  highly  desirable  for  achieving  better  control  over  the  far-held  profile  by  minimizing  the 
blurring  effects  due  to  phase  aberration. 

The  structure  can  also  be  reduced  by  using  optical  relaying,  which  can  be  readily  imple¬ 
mented  by  the  configuration  shown  in  Fig.  10.  Here  one  replaces  the  diverging  collimator 
of  Fig.  2  by  a  converging  mirror  located  a  distance  equal  to  its  focal  length  fc  >  0  be¬ 
yond  the  beam  focus,  thus  giving  an  aperture-collimator  separation  distance  Zac  =  f a  +  fc  ■ 
The  scaled  distance  now  reverses  sign  to  become  Zac  —■ *  —Xzac/  (DaDc)  ,  but  because 
its  magnitude  is  the  same  as  in  the  previous  configuration,  the  structure  at  the  collima¬ 
tor  remains  qualitatively  similar  to  that  of  Fig.  6a.  In  this  configuration,  an  image  I\  of 
the  amplifier  aperture  forms  at  a  distance  zci  =  zacDc  /  Da  =  13.25  m  beyond  the  colli¬ 
mator.  (This  follows  from  the  requirement  that  the  scaled  propagation  distance  vanishes: 
Zai  =  Zac  +  Zci  =  —Xzac/  ( DaDc )  +  Xzci/D q  =  0.)  At  plane  R\  located  a  distance 
Z(j  /  beyond  A,  the  scaled  propagation  distance  increases  to  Z(:  r  =  Zac  I ,  so  the  transverse 
fluence  distribution  becomes  an  exact  inverted  image  of  Fig.  6a.  The  amplitude  of  the 
structure  begins  to  increase  beyond  R\,  so  that  would  be  the  last  chance  to  place  the  first 
optic  of  a  1  :  1  image  relay  telescope  without  subjecting  it  to  larger  fluence  spikes.  The 
structure  at  the  second  optic  R->  (at  a  distance  2zci  from  R \ )  is  comparable  to  that  of  R\ . 
and  the  telescope  creates  a  second  image  plane  R  at  a  distance  4 zci  =  53  m  from  the  first 
one.  Because  the  final  focusing  lens  L  could  be  safely  placed  an  additional  distance  of  only 
zci  beyond  this  image,  the  largest  total  available  delay  time  beyond  I\  is  5 zci/c  =  221  ns. 
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This  would  be  insufficient  for  some  of  the  beams  in  a  large  angularly-multiplexed  system,  so 
these  beams  would  require  an  additional  relay  telescope. 

In  principle,  it  would  also  be  possible  to  reduce  the  structure  by  taking  advantage  of 
the  ISI  filtering  process  discussed  above.  In  an  angularly-multiplexed  system,  however, 
that  approach  would  require  long  propagation  paths  to  the  time  demultiplexing  mirrors,  as 
well  as  to  the  focusing  lenses.  It  would  result  in  excessively  long  propagation  paths  for  all 
the  beams,  which  could,  in  turn,  make  beam  alignment  more  difficult  and  create  possible 
imprinting  problems  at  the  far-held. 

Summary  and  Conclusion 

Although  ISI  provides  ultrasmooth  illumination  at  the  far-held  of  a  laser,  it  still  allows 
speckle  to  develop  in  the  quasi  near-held  beams.  This  paper  has  compared  simulations  of 
that  speckle  with  measurements  in  the  output  beams  of  the  Nike  laser  and  explored  options 
for  controlling  the  problem  in  future  KrF  drivers.  The  speckle  arises  from  phase  aberration 
in  the  system,  and  builds  up  as  the  beams  propagate  away  from  the  exit  pupil  plane  at  the 
hnal  amplifier.  Because  of  the  spatial  incoherence  of  ISI  light,  the  speckle  is  significantly 
smaller  than  that  developed  by  coherent  light;  nevertheless,  it  remains  a  damage  issue  for  the 
higher  huences  and  long  delay  paths  required  in  future  angularly-multiplexed  KrF  drivers. 

We  have  examined  two  techniques  for  calculating  the  average  intensities.  The  first  prop¬ 
agates  a  randomly-selected  stochastic  optical  held  amplitude  through  the  laser  system  to 
evaluate  the  instantaneous  intensity,  then  repeats  this  procedure  with  multiple  new  statisti¬ 
cal  realizations  to  accumulate  the  time-average  intensity.  The  second,  which  is  much  faster, 
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evaluates  the  ensemble-average  intensity  envelope  directly  using  FFT  and  a  correlation  func¬ 
tion  algorithm.  Benchmark  tests  show  excellent  agreement  (within  <  1%  rms)  between 
the  two  approaches,  providing  conhdence  that  the  algorithm  can  be  routinely  applied  to 
calculations  and  analysis. 

The  Nike  simulations,  which  are  in  good  qualitative  agreement  with  the  photochromic 
him  measurements,  show  that  as  the  distance  from  the  amplifier  increases,  (a)  the  structure 
becomes  coarser  and  (b)  the  hot-spot  amplitudes  increase  at  first,  then  appear  to  saturate 
and  slowly  decrease.  We  also  performed  simulations  to  explore  options  for  reducing  the 
near-held  structure  in  future  KrF  drivers.  One  possibility  is  to  reduce  the  phase  aberration; 
e.g.,  a  2 -fold  reduction  in  point-spread  width  would  also  reduce  the  amplitudes  by  ~  2- fold. 
The  second  option  uses  optical  relaying,  but  that  would  probably  be  expensive  to  implement 
in  some  of  the  longer  beam  paths  required  for  demultiplexing. 

Acknowledgement:  This  work  is  supported  by  the  USDOE 

Appendix 

This  appendix  examines  the  condition  necessary  to  satisfy  the  quasi-monochromatic  ap¬ 
proximation.  The  basic  criterion  requires  the  optical  bandwidth  Au  to  be  narrow  enough 
that  it  does  not  suppress  interference  structure  in  the  propagation.  To  get  interference  at 
point  between  light  originating  from  the  two  aperture  points  and  X'A  (e.g.,  see  Fig. 
2),  the  total  path  difference  between  the  two  rays  must  be  smaller  than  the  coherence  length 
crc  ~  cj  Av  ~  300  /iin  for  Au  ~  1  THz.  In  the  paraxial  approximation,  the  pathlength  for  a 
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single  ray  is  the  total  phase  shift  divided  by  k0  =  2ir/\: 

l  (XA,  XL)  ~  ^AC  +  %cl  +  (27T/A)-1  (tt /Zal)  \Xl  -  XA|2  ,  (Al) 

where  the  transverse  contribution  (7 t/Zal)  |X^  —  XA|2  comes  from  expression  (7).  (As  ex¬ 
pected,  the  minimum  value  of  l  (XA,X^)  occurs  for  rays  that  satisfy  the  geometrical  optics 
condition  X^  =  XA.)  The  criterion  1 1  (XA,X^)  —  l  (XA,Xx,)|  <  c/Av  then  becomes 

Z~AlL  |  |XL  -  XA|2  -  |XL  -  X'a|2|  <  2z/0/ Av ,  (A2) 

where  v0  =  c/A  ~  1200  THz.  For  spatially- coherent  light  the  largest  factor  in  (A2)  would 
occur,  for  example,  when  XA  =  — XA  =  X^  and  |X^|  =  l/\/2;  the  criterion  then  reduces  to 
Zal  >  Av/vq  ~  8xl0~4.  In  Figure  8,  this  is  marginally  satished  at  the  target  chamber,  where 
ZAl  =  ZAc+Zcl  =  9.84xl0~4,  but  not  at  the  recollimator,  where  Zal  —■ *  Zac  =  1.46xl0~4. 
Thus  the  very  high  spatial  frequency  ringing  shown  at  the  recollimator  in  Fig.  8  (i.e.,  the  fine 
Fresnel  structure  at  the  edges)  would  actually  be  washed  out  at  the  bandwidths  of  interest. 
For  ISI  light  of  far-held  XDL  width  Np,  the  factors  in  (A2)  would  be  limited  by  the  diameter 
1/Np  <C  1  of  the  spatial  coherence  zones;  hence  the  criterion  becomes 

Zal  >  (2A^)  2  Azz/z/0;  (A3) 

which  is  easily  satished  even  at  the  recollimators  when  NF  =  75  XDL. 
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Figure  Captions 


Fig.  1.  Nike  amplifier  staging  and  angularly  multiplexed  optical  system.  The  object  aper¬ 
ture  is  uniformly  illuminated  by  spatially  and  temporally  incoherent  light,  then  imaged 
onto  the  target  through  the  amplifier  system.  The  beams  overlap  at  the  amplifiers, 
which  are  placed  at  or  near  images  of  the  object’s  Fourier-transform  (pupil)  plane. 

Fig.  2.  Schematic  of  one  beam  of  the  Nike  optical  system  starting  at  the  60  cm  amplifier 
(aperture  stop  A),  converging  to  the  recollimating  mirror  C,  then  propagating  along 
delay  path  zcl  =  %cm  +  zml  to  focusing  lens  L.  The  wavy  line  at  A  represents  the 
random  phase  of  complex  transmission  TA(x).  The  amplifier  is  double-passed,  so  the 
large  lens  represents  a  concave  mirror.  The  dotted  lines  near  C  and  L  illustrate  the 
placement  of  the  film  packs  where  the  beam  #10  profiles  were  measured. 

Fig.  3.  Simulations  of  the  point-spread  function  and  ISI  focal  profile  on  the  Nike  laser, 
assuming  a  75  XDL  wide  ”  top-hat”  object  profile,  a  uniformly  illuminated  60  cm 
amplifier,  and  random-phase  aberration  of  A/4  rms  amplitude  and  ~  15  XDL  point- 
spread  width. 

Fig.  4.  Comparison  of  ISI  far-field  and  near-field  intensities  averaged  over  8000  coherence 
times  (thin  lines)  with  envelope  intensities  (heavy  lines).  In  these  simulations,  the 
object  fluence  is  a  32  XDL  wide  ” top-hat”  profile  and  the  random-phase  aberration  is 
modelled  by  a  ~  7  XDL  point-spread  width. 
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Fig.  5.  Photochromatic  film  images  and  calibrated  lineouts  at  the  Nike  laser  recollimator 
mirror  array  and  target  chamber  windows:  (a)  Beam  #10  at  the  recollimator  mirror 
and  (b),  (c)  Beams  #38  and  #10,  respectively,  at  the  target  chamber. 

Fig.  6.  Simulations  of  ISI  intensity  envelope  images  and  lineouts  at  the  Nike  laser  recolli- 
mator  array  and  target  chamber  windows:  (a)  Beam  #10  (or  #38)  at  the  recollimator 
mirror  and  (b),  (c)  Beams  #38  and  #10,  respectively,  at  the  target  chamber. 

Fig.  7.  ID  mode  spectra  of  measured  ( - )  and  simulated  ( - )  near- field  fluence 

structure  at  the  Nike  laser  recollimator  array  and  target  chamber  windows:  (a)  Beam 
#10  recollimator  mirror  and  (b),  (c)  Beams  #38  and  #10,  respectively,  at  the  target 
chamber.  The  measured  data  has  larger  statistical  variations  because  it  averages  over 
fewer  horizontal  lines. 

Fig.  8.  Simulations  of  coherent  intensity  images  and  lineouts  at  the  recollimator  mirror 
and  target  chamber  window  of  Nike  beam  #10.  (Note  the  extended  ordinate  scale  in 
the  lineouts,  compared  to  the  earlier  figures.) 

Fig.  9.  Same  as  Fig.  6a, b,  except  that  the  point-spread  width  is  8  XDL  instead  of  15  XDL 

Fig.  10.  Schematic  of  an  alternate  ISI  optical  system  at  and  beyond  the  final  aperture 
stop  A,  showing  the  beam  focusing  before  the  recollimating  mirror  C,  propagating 
through  the  first  amplifier  image  plane  I\  to  a  1  :  1  relay  telescope  beginning  at  Ri, 
and  producing  a  second  amplifier  image  plane  J2  before  the  final  focusing  lens. 
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Fig.  1:  Nike  amplifier  staging  and  angularly  multiplexed  optical  system.  The  object  aperture  is  uniformly  illuminated 
by  spatially  and  temporally  incoherent  light,  then  imaged  onto  the  target  through  the  amplifier  system.  The  beams 
overlap  at  the  amplifiers,  which  are  placed  at  or  near  images  of  the  object’s  Fourier- transform  (pupil)  plane. 
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Fig.  2:  Schematic  of  one  beam  of  the  Nike  optical  system  starting  at  the  60  cm  amplifier  (aperture  stop  A),  converging  to  the 
recollimating  mirror  C,  then  propagating  along  delay  path  zCL  =  zCM  +  zML  to  focusing  lens  L.  The  wavy  line  at  A  represents  the 
random  phase  of  complex  transmission  T4(x).  The  amplifier  is  double-passed,  so  the  large  lens  represents  a  concave  mirror.  The 
dotted  lines  near  C  and  L  illustrate  the  placement  of  the  film  packs  where  the  beam  profiles  were  measured. 
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Fig.  3.  Simulations  of  ISI  on  the  Nike  laser  with  a  75  XDL  wide  “top  hat”  object  profile,  uniformly  illuminated  60  cm  amplifier, 
and  random-phase  aberration  of  A/4  rms  amplitude  giving  a  point-spread  function  of -15  XDL  width. 
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Fig.  4.  Comparison  of  ISI  intensities  averaged  over  8000  coherence  times  (thin  lines)  with  envelope  intensities 
(heavy  lines)  at  the  far- (i eld,  the  laser  aperture,  and  final  focusing  lens.  In  these  simulations,  the  object  fluence  is  a 
32  XDL  wide  “top-hat”  profile  and  the  random-phase  aberration  is  modelled  by  a  ~7  XDL  point-spread  width. 
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Fig.  5.  Photochromatic  film  images  and  calibrated  lineouts  at  the  Nike  laser  recollimator  mirror  array  and  target  chamber 
windows:  (a)  Beam  #10  at  the  recollimator  mirror  and  (b),  (c)  Beams  #38  and  #10,  respectively,  at  the  target  chamber. 
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Fig.  6.  Simulations  of  ISI  intensity  envelope  images  and  lineouts  at  the  Nike  laser  recollimator  array  and  target  chamber 
windows:  (a)  Beam  #10  (or  #38)  at  the  recollimator  mirror  and  (b),  (c)  Beams  #38  and  #10,  respectively,  at  the  target  chamber. 
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Fig.  7.  ID  mode  spectra  of  measured  ( - )  and  simulated  ( - )  near-field  fluence  structure  at  the  Nike  laser  recollimator 

array  and  target  chamber  windows:  (a)  Beam  #10  recollimator  mirror  and  (b),  (c)  Beams  #38  and  #10,  respectively,  at  the 
target  chamber.  The  measured  data  has  a  larger  statistical  variations  because  it  averages  over  fewer  horizontal  lines. 
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Fig.  8.  Simulations  of  coherent  intensity  images  and  lineouts  at  the  recollimator  mirror  and  target  chamber 
window  of  Nike  beam  #10.  (Note  the  extended  ordinate  scale  in  the  lineouts,  compared  to  the  earlier  figures.) 
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Fig.  9.  Same  as  Fig.  6a, b,  except  that  the  point-spread  width  is  8  XDL  instead  of  15  XDL 


Fig.  10.  Schematic  of  an  alternate  ISI  optical  system  at  and  beyond  the  final  aperture  stop  A,  showing  the  beam  focusing  before 
the  recollimating  mirror  C,  propagating  through  the  first  amplifier  image  plane  /,  to  a  l  :  1  relay  telescope  beginning  at  Rh  and 
producing  a  second  amplifier  image  plane  I2  before  the  final  focusing  lens. 


